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We consider a Friedmann-Robertson-Walker cosmological model with some exotic perfect fluid 
with spin known as the Weyssenhoff fluid. The possibility that the dark energy may be described 
in part by the Weyssenhoff fluid is discussed. The observational constraint coming from supernovae 
type la observations is established. This result indicates that, whereas the cosmological constant 
is still needed to explain current observations, the model with spin fluid is admissible. For high 
redshifts z > f the differences between the model with spin fluid and the cold dark matter model 
with a cosmological constant become detectable observationally for the flat case with £lm,o = 0.3. 
From the maximum likelihood method we obtain the value of f2 s ,o = 0.004 ± O.Of 6. This gives us 
the limit Q St o > —O.Of 2 at the la level. While the model with "brane effects" is preferred by the 
supernovae la data, the model with spin fluid is statistically admissible. For comparison, the limit 
on the spin fluid coming from cosmic microwave background anisotropies is also obtained. The 
uncertainties in the location of a first peak give the interval -1.4 x f0~ 10 < O Sj0 < -10~ 10 . From 
big bang nucleosynthesis we obtain the strongest limit f2 s ,o ^ — 10 -20 . The interconnection between 
the model considered and brane models is also pointed out. 



I. INTRODUCTION 

In 1923 Cartan introduced the intrinsic angular momentum in the theory of relativity (as a classical quantity) 
before it was introduced as spin into quantum theory by Goudsmit and Uhlenbeck in 1925. The classical spin can be 
introduced in general relativity in two distinct ways. The first one is to introduce spin as a dynamical quantity in 
special and then in general relativity without changing the geometry, i.e., without modifying the metric of spacetime 
0. [1 [1113 El The spin introduced in this way showed more or less a similarity to the spin of quantum mechanics (and 
the Dirac theory of the electron). The second, more satisfactory way to introduce the intrinsic angular momentum 
is by generalizing the structure of spacetime. It was done by Cartan by assuming the metric and the non-symmetric 
affine connection as independent quantities and was further developed by Hehl 7], Trautman Q and Kopczynski 
0, IToj . This assumption allowed a definition of the torsion of spacetime and its connection with the torsion with 
spin. In the framework of Einstein-Cartan theory for the Friedmann universe Trautman |ll| made the conclusion that 
torsion avoids the singularity and stops the collapse in closed models at the moment of minimum radius about 1 cm 
with matter density p ~ 10 55 g cm~ 3 
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Let us note that the above formulas are valid for chaotic spin distribution na. The effects of spin fluid are important 
in a low-energy limit of the superstring theory (see for general discussion |13|) which is supergravity whose integral 
part is torsion. 

The torsion contributes to the energy-momentum of spin fluid which has the form [T^| 

T jS = u n u AP + P- 2s2 ) - 9^{P - s 2 ) 

with s 2 = s^vS^" , where the spin s leads to an effective negative pressure and eliminates the singularity. 

Let us consider a world model with the Robertson- Walker symmetries which is filled with 'perfect fluid with spin'. 
As it is well known [l5| the macroscopic spin tensor r x may be expressed in terms of the spin density tensor and 
four- velocity of the fluid u M (u^u^ = — 1) 

t£, = u x S^ S^u" = 0. (1) 

To describe the material content of the considered model we use the hydrodynamical description in terms of the 
energy-momentum tensor which in the general relativity limit reduces to perfect fluid characterized by the energy 
density e and the isotropic pressure p. In analogy, the physical content of the model in the Einstein-Cartan theory, 
based on the classical description of spin with equation Q , may be called "perfect fluid with spin" . This is surely 
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the simplest type of hydrodynamic continuum of use for our aims. It is an extension of the well known semi-classical 
model of spin fluid from special relativity; we call it the Weyssenhoff fluid [|| . The influence of the macroscopic spin 
present in the fluid on the dynamics of the universe is then described by contributions to the energy density and 
pressure 

e c f{ = e-^S 2 , p cS = p-^S 2 , (2) 

where S 2 = ^S^S^ , and S^ v is th e sp in density tenso r Il6| . 

Supernovae type la observations |l7l ITU Il9l I20I l2ll l22j indicated that the Universe's expansion has started to 
accelerate during recent cosmological times. These observations (as well as cosmic microwave background observations) 
suggest that the energy density of the Universe is dominated by a dark energy component with negative pressure, 
driving the acceleration. The most natural candidate to represent dark energy is the cosmological constant. However, 
it is necessary to require a fine tuning of 120 orders of magnitude in order to obtain agreement with observations |23| . 

In this work, we discuss the possibility that some part of the dark energy may be represented by the perfect fluid 
with spin, which is characterized by an equation of state in form Ij2fl. 

Let us consider the dust fluid (p = 0) of particles with spin 1/2 and mass m. The energy density e and absolute 
value of the spin density S depend on number n of particles in a unit volume 

e = nm, b = —n. 

Hence we have 

Therefore, effective pressure is negative. Because the number of particles in comoving volume n(t) oc a~ 3 (i), where 
a is the scale factor and e = e mj oa (t) is the energy density of dust, the energy density e and pressure p are never 
greater than 
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respectively. The energy density of the spin is negative but e e g is assumed to be positive. The ratio of the spin-induced 
term to the standard energy term is 

(8irGS) 2 _ 8wGh 2 p 
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With nucleons as dust particles, the above ratio is of the order 8 x 10~ 55 p (where p is expressed in the cgs units). 
Although the spin-spin contact interaction term appears to be negligibly small even in what may be regarded now as 
superdense matter, it may play an essential role when the ratio approaches unity, i.e., in the earliest stages of evolution 
of the Universe. The spin-spin term produces something which may be called after Kopczynski the "centrifugal force" 
and which is able to prevent the occurrence of singularities in cosmology 0, ITol ] . 

Let us note now that effects of spin are dynamically equivalent to introducing into the model some additional 
non-interacting fluid for which the equation of state is 

p s = w s p s , (4) 

where w s = 1, p s = p s .oa~ 6 , p Sj o = — (fi 2 /16)n(0), just as for the Zeldovich stiff matter (or brane effects with dust on 
a brane with negative tension). 



II. THE FRW MODEL WITH SPIN ON TWO-DIMENSIONAL PHASE PLANE 

The dynamics of the Friedmann-Roberston- Walker (FRW) model with the Weyssenhoff fluid can be represented on 
the phase plane (H, e), where H = d(lna)/dt is Hubble's function, in the following way 

-tt = ~H 2 - - (e eH + 3 PcS ) + - 5a 
dt 6 3 
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FIG. 1: The phase portrait of system ©. There are the four critical points: stable node, unstable node, saddle point and 
center. The presence of the center in the finite region indicates that the model is structurally unstable. 



where p c s and e c g are parameterized by the energy density of dust (see formula ©). Of course, system (|5j) has the 
first integral 

£cff -3i/ 2 = ^-A, (6) 

where k G {0, ±1} is the curvature constant. Additionally, we have the conservation condition for non-interacting 
spin fluid with energy density e 

f = m 

The critical points of (J5J can be one of two admissible types: static if Ho = 0, eo ^ 0, or non-static if Hq ^ and 
eo = 0. In the first case critical points lie on the intersection H-axis with the boundary condition e e g + 3p c ff = (the 
cosmological constant is formally included into p e g and /o e ff in the standard way p\ = A and p\ = —A). The second 
type of critical points lie on the intersection of e-axis with the trajectory for flat model {e = 3H 2 }. 
The physically admissible domain for trajectories is 

{(H, e) : ffeM and < e < e max , e max : e cff (e max ) = -A}. 

The phase portrait of system 10) with A > is demonstrated on Fig. ^ 

The differences in the behavior of trajectories are manifest at high densities. Then the structure of dynamical 
behavior at infinity is modified. To illustrate the behavior of trajectories at infinity, system (JSJ) is represented on 
Fig. m the projective coordinates 

z = Jt> u ^~h_' (z,u)-map (8a) 

I H 

v = -, w = — , (v,w) — map. (8b) 



The two maps (z, u): z — Q, —oo < u < oo and (u, w) : v = 0, — oo < w < oo cover the behavior of trajectories at the 
infinity circles H = oo and e = oo. 




FIG. 2: The phase portraits of system © in the projective coordinates (z,u) and (v, w). From the portrait in the former 
coordinates we can see that there is a saddle point in the infinity. 
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III. DYNAMICS OF THE MODEL WITH SPIN FLUID 



It is well known that the first integral of the FRW equation can be used to construct a Hamiltonian function. We 
take advantage of this feature in the considered model. The integration of equation (|5b|) gives p — p$a ~ 3 . Therefore 
the right-hand side of the Raychaudhuri equation (|5al) can be expressed in terms of the scale factor a(t) as 
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Equation 10 can be rewritten in a form analogous to the Newton equation of motion in the one-dimensional config- 
uration space {a: a 6 K+} 



where the potential function 



dV(a) 
da 



V(a) = -le^aa- 1 + [ 'i ) - ^Aa 2 + V 
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and Vq — const. First integral © should correspond with the Hamiltonian integral of motion, hence 
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(10) 
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V(a) + Y = V --. 



Now we construct the Hamiltonian function 
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and then the trajectories of the system lie on the energy level Ti = E, but if we choose Vq = k/2 then the physical 
trajectories lie on the zero-energy level H = E = 0, which coincides with the form of first integral 10. Then finally 



we obtain 



or in general 



where 
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and n(Q) is an initial number of particles in the unit comoving volume. 

Finally we obtain dynamics reduced to the form of particle like problem in the one-dimensional potential 



x = y 



with 



dV 

dx 



V(x) = --(e mfl x 1 +e ai ox 4 - e A fiX 2 - 3k) 



where x — a, y = d and the above system has the first integral 



(13) 
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(15a) 
(15b) 

(16) 



+ V(x) = 0. 



(17) 
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In order, basic dynamical system (|15|) is then rewritten as 

Y = ^,° + lJ2 Q *>° v ~ {1+3Wi) ( 18a ) 

i 

v = -^Y,n ii0 {l + 3w i )v- ( - 2+3w ^ (18b) 

i 

where v = a/do, t — ► T = |i?o|i is a new time variable denoted as dot in (|18fl : € {^ m ,0j ^s,o> ^A,o} are density 
parameters, fii = Pi/ p cri tyo — Pi.o/iH^ and the subscript means that a quantity with this subscript is evaluated 
today (at time to); pi — Pi,oa~^ 1+Wi \ where w m = (dust), w s — 1 (spin fluid), w\ = —1 (cosmological constant), 
and pi =Wipi. 

The representation of dynamics as a one-dimensional Hamiltonian flow allows to make the classification of possible 
evolution paths in the configuration space which is complementary to phase diagrams. It also makes simpler to discuss 
the physical content of the model. Finally, the construction of the Hamiltonian allows to study quantum cosmology 
models with spin fluid in full analogy to what is usually done in general relativity. 

From equation l|14f> we can observe that trajectories are integrable in quadratures. Namely, from the Hamiltonian 
constraint Ti = E = we obtain the integral 

r da 

t-t = / - 19 

j ao v^sv^y 

For some specific forms of potential function l|16(l we can obtain exact solutions. 

It is possible to make the classification of qualitative evolution paths by analyzing the characteristic curve which 
represents the boundary equation of the domain admissible for motion. For this purpose we consider the equation of 
zero velocity, a — which constitutes the boundary M = {a: V(a) = 0}. 

From equation l|16(l the cosmological constant can be expressed as a function of x as follows 

A(x) = ^(e^oaT 1 - e S:0 x~ 4 - 3fc). (20) 

The plot of A(x) for different k is shown in Fig.0 The domain under the curve A(cc) is non-physical, and we consider 
the evolution path as a level of A = const and we classify all evolution models with respect to their quantitative 
properties of dynamics. For negative A there are oscillating solutions without a singularity. 

The next advantage of representing dynamics in terms of Hamiltonian is possibility to discuss how trajectories 
along which the acceleration condition a — —dV/da > is satisfied arc distributed in the phase space. One can easily 
observe this phenomenon from the geometry of the potential function. In the phase plane the area of acceleration is 
determined by y > or by the condition that V(x) is a decreasing function of its parameter 

-7TOn,0^ 2 -2ft s ,0V~ 5 + n Afi V>0 (21) 



where 



— , a,o < 0, M = o. 
a 



Independent observations of supernovae type la made by the Supernovae Cosmology Project and High-z Survey 
Team indicate that our Universe is presently accelerating. There is a fundamental problem in explaining this accel- 
eration. If we introduce the cosmological constant and assume ftk,o = (the universe is flat) then the best fit model 
is fl Afi = 0.72, Sl m>0 = 0.28. 

The effects of spin fluid cannot dominate the matter contributions during the whole evolution of the universe. But 
we argue that the effects of spin degree of freedom are important in early universe. In any case they should be smaller 
than or comparable with the matter contribution because 

^ = tt m ,o^- 3(1+7 + a.o^- 6 > o. 

Ho 

The formalism presented gives us a natural base to discuss this problem for the FRW model with the Weyssenhoff 
fluid. It is convenient to introduce a new variable z = v~ 3 , then inequality l|21|l reduces to the quadratic inequality 

-^mfiZ + 2|fL 5i0 |z 2 + Ov i0 > 0. 
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FIG. 3: The dependence of A on a. 



Therefore the Universe is presently accelerating provided that 

-iftm,0 + 2|fi s , | +^A,0 > 0. 

For instance, for fl m .o = 0.3, Qa,o = 0, |fi s ,o| > j^m.o is required. Therefore, present experimental estimates based 
on baryons in clusters on one hand giving f2 m ,o ~ 0.3 and the location of the first peak in the cosmological microwave 
background detected by Boomerang and Maxima suggesting an early, filled by spin matter, flat universe on the other 
hand, imply that our Universe without a cosmological term would be presently accelerating only if 




The required value of |f2 s ,o| seems to be unrealistic (see the next section) and therefore the cosmological constant 
term is still needed to explain the present acceleration of the Universe. As we will see in the subsequent analysis the 
effect of spin fluid is negligible in the present epoch. Therefore the spin is an additional factor that can influence the 
dynamics of the early universe together with for example the cosmological constant. 



IV. REDSHIFT-MAGNITUDE RELATION FOR THE MODEL WITH SPIN FLUID 

Cosmic distance measures like the luminosity distance depend sensitively on the spatial geometry (curvature) and 
the dynamics. Therefore, the luminosity depends on the present density parameters of different components and their 
form of the equation of state. For this reason the redshift-magnitude relation for distant galaxies is proposed as a 
potential test for the FRW model with spin fluid. 

Let us consider an observer located at r — at the moment t — t which receives a light ray emitted at t = t\ from 
the source of the absolute luminosity L located at the radial distance n . The redshift z of the source is related to the 
scale factor at the two moments of evolution by 1 + z = a(to)/a(ti) — ao/a. If the apparent luminosity of the source 
as measured by the observer is I then the luminosity distance 4 of the source is defined by the relation 
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where cLl = (1 + z)aor\. 

The important test to verify whether the spin fluid FRW model may represent dark energy is the comparison with 
the supernova type la data. In order to do so, we evaluate the luminosity distance in the flat FRW model with spin 
fluid. In such a model, the luminosity distance reads 



d L {z) 



l + z 



dz' 
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From this expression the following relation between the apparent magnitude m and absolute magnitude M is obtained 

dz' 



m — M = 5 log 



(l + z) 



o y/n m ,o{i + z'f + n s . (i + z'f + n Afi 



where M = M - 5 log H + 25. 

In order to compare with the supernova data, we compute the distance modulus 

^ = 51og(d L ) + 25 

where dh is in Mps. The goodness of fit is characterized by the parameter 
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°/*0,» + a fiz,i 



(22) 



In expression (|2"2"|) . is the measured value, is the value calculated in the model described above, a^ oi is the 
measurement error, cr 2 i is the dispersion in the distance modulus due to peculiar velocities of galaxies. 



V. THE MODEL WITH SPIN FLUID TESTED BY SUPERNOVAE 



We test the model with spin fluid using sample A of the Perlmutter SN la data. In order to avoid any possible 
selection effects we work with this full sample of 60 supernovae. In the statistical analysis we use the maximum 
likelihood method with the marginalization procedure |l9j . 

First we estimated M. from the full sample of 60 supernovae. For the flat model we obtained M. = —3.39. 

The result of the statistical analysis is presented in the figures. Figure 01 illustrates the confidence level as a function 
of £l m Q, f2 Sj o for the flat model (Hk,Q = 0) minimized over M with £Ia,o — 1 — ^m,o ™ ^fe,o ~ ^s,o- 

In the case of Qk,o 7^ 0, the preferred values of the pairs (r2 m ,0j ^A,o) are shown in a standard way after minimizing 
over M. on Fig. One can see that the non-zero cosmological constant is still required. 

Applying the marginalization procedure over Qk,o = (~l> 1)> ^s,o < 0, A4 G (—4, —3) we find the lowest value of % 2 
for each pair of values (£l m ,o, ^A,o) as shown on Fig. This figure shows the favored region of values of (f2 m ,0j ^A.o) 
(best fit). 

Figure [3 shows the density distribution for fi Sj o in the flat model with flmfl — 0.3. This distribution is obtained 
from the marginalization over M. . The positive values of f2 Si o can be interpreted as the brane effects with dust on the 
brane [24^. One can see that Q s ,o > —0.012 at the confidence level la and f2 Sj o > —0.026 at the confidence level 2a. 
The limit value of il s o — —0.012 is used in our further analysis. 

In Fig. [S] we present the plot of residuals of redshift-magnitude relationship for the supernovae data. With the 
increasing impact of spin (lower f2 s ,o) the high-redshift supernovae should be fainter than the expected by the ACDM 
model. For Slg.o = —0.012 the difference between the spin model and the ACDM model should be detectable for z > 1. 

The angular diameter of a galaxy is defined as 

d(z + lf 

where d is the linear size of a galaxy. In the standard cosmology the flat dust-filled universe 9 has the minimum 
value z m in = 5/4. From Fig. [5] we can see that, if the spin matter is present, its influence on the predicted 9 is weak. 
Theoretically it is possible to test values of r2 s ,o from the angular diameter minimum value test. But because of 
evolutionary effects and observational difficulties the predicted differences are too small to be detected. 
Now let us briefly discuss the effect of spin fluid on the age of the Universe which is given by 

^ _ 1 ^ 2f2 m ,o + 2Qa,o + ^/^a,o 
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FIG. 4: Confidence levels on the plane (fi m ,o, f2 s ,o) minimized over M for the flat model, and with J7a,o = 1 — f2m,o — Ofe.o — ^s,o • 
The figure shows the ellipses of the preferred value of f2 m ,o and Qa,o- The results prefer the positive value of S7 s ,o, while the 
negative values are allowed (i.e., spin fluid can exist). 
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FIG. 5: Confidence levels on the plane (fi m ,o,^A,o) minimized over M, and with £1a,o = 1 — ^m,o — Qk,o — ^s,o- The figure 
shows the ellipses of the preferred values of f2 m ,o and Qa,o- 



10 




FIG. 6: Levels of constant \ 2 on the plane (fi m ,o,fiA,o) marginalized over M, and with S7a,o = 1 — f2 m ,o — f2fe,o — ^a,o The 
figure shows the preferred value of fi m ,o and S1a,o- 

In Fig. we plot the age of the Universe in Gyr for the flat model for different Q s ,o- Taking f2 m) o = 0.3 we obtain 
for Sl s ,o = —0.013: to = 13.347, and for Q s ,o = 0: to = 14.436. We can see that the spin fluid lowers significantly the 
age of the Universe. The cosmological constant is still needed to explain the problem of the age of the Universe. 



VI. CMB PEAKS IN THE MODEL WITH SPIN FLUID 



The cosmic microwave background (CMB) peaks arise from acoustic oscillations of the primeval plasma. Physically 
these oscillations represent hot and cold spots. Thus, a wave that has a density maximum at the time of last scattering 
corresponds to a peak in the power spectrum. In the Legendre multipole space this corresponds to the angle subtended 
by the sound horizon at the last scattering. Higher harmonics of the principal oscillations which have oscillated more 
than once correspond to secondary peaks. 

For our end it is very important that the locations of these peaks are very sensitive to the variations in the model 
parameters. Therefore, it can serve as a sensitive probe to constrain the cosmological parameters and discriminate 
among various models. 

The locations of the peaks are set by the acoustic scale I a which can be defined as the angle 9 a subtended by the 
sound horizon at the last scattering surface. The acoustic scale I a = k/Qa in the flat model is given by 

f^dec dz 

j JO H(z') 

La - 77 r°° ' dz~ 
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FIG. 7: The density distribution for Q Bl o in the model with spin effects (Q s ,o < 0) and brane effects (fi s ,o > 0). We obtain that 
fis.o > —0.012 at the confidence level la and fi s ,o > —0.026 on the confidence level 2a. For simplicity, we only mark limits for 
the spin side. 



where 



H(z) = ff ^a-,o(i + zf + n s , (i + zf + n m , (i + zf + ft A , 



and c s is the speed of sound in the plasma (we assume additionally the presence of radiation in the model). The 
sound velocity can be calculated from the formula 



dp c 



_ 6n B , (l + z) 3 + fn 7 , (l + z) 

Us ~ dp cfi ~ 3fJ m , + 4fi 7 , (l + z) + 6fJ S(0 (l + zf ' 

In the model of primeval plasma, there is a simple relation 

lm ~lA{m- <j) m ) 



between the location of m-th peak and the acoustic scale J25, 26]. The prior assumptions in our calculations are as 
follows n r ,o = 9.89 x 10~ 5 , f2b,o = 0.05, and the spectral index for initial density perturbations n=l, and h = 0.65. 

The phase shift is caused by the prerecombination physics (plasma driving effect) and, hence, has no significant 
contribution from the term containing spin in that epoch. Because of above assumptions the phase shift cj> m can be 
taken from standard cosmology 25] 
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FIG. 8: For f2 m ,o = 0.3, residuals between the Einstein-de Sitter model and three cases: the Einstein-de Sitter itself (zero line), 
the ACDM flat model (middle curve), and the flat model with spin matter for f2 s ,o = —0.012 (highest curve). 

where il b ,oh 2 = 0.02, r(z doc ) = p T (zd ec ) / Pm(zdec) = fi r ,o(l + ^doc)/^m,o, ^r,o = ^7,0 + ^,0, ^7,0 = 2.48/i~ 2 x 10~ 5 , 
= l-7/i~ 2 x 10~ 5 , r(zdcc) is the ratio of radiation to matter densities at the surface of last scattering. 
The influence of spin on the location of the peaks is to shift them towards higher values of I. For example, for 
^m,o = 0.3, f2b,o = 0.05, h — 0.65, the different choices of f2 Sj o yield the following 

^s,0 = 10 10 : ^peak.l = 186, 'peak,2 

^s,0 = 1-4 X 10 10 : £peak,l = 223, ^peak,2 

On the other hand from the Boomerang observations [27j we obtain i pea k,l = 200 — 223, / pC ak.2 
uncertainties in values l pe ak can be used in constraining cosmology with spin fluid, namely 

10~ 10 < |fi s , | < 1.4 x 10~ 10 

from the location of the first peak. 

VII. BBN IN THE MODEL WITH SPIN FLUID 

We consider the non-relativistic matter with density p. In this case the spin effects scale like —(1 + z) 6 . It is clear 
that such a term can lead to accelerated expansion and the detailed analysis of supernovae la data requires it to be 
very subdominant today. 

However, going back in time to the big bang nucleosynthesis (BBN) period, the term fi Si o = —0.012 would be 
dominant at redshift z > 2. In such a model, radiation domination would never occur and all the BBN predictions 
would be lost. In practice, the preferred value obtained from SN la data gives the term p 2 which is far too large to 
be compatible with the BBN which is very well tested area of cosmology and does not allow for significant deviation 
from the standard expansion law, apart from very early times before the onset of the BBN. The consistency with the 
BBN seems to be crucial issue in spin fluid cosmology j23, IS] . For this reason, not to suffer from the contradiction 
with the BBN, the contribution of spin fluid f2 Si0 cannot dominate over the standard radiation term before the onset 
of BBN, i.e., for z = 10 8 



= 441 
= 530 

= 509-561. Therefore, 



-f7 Sj0 (l + z) 6 < f7 r ,o(l + z) 4 and -O s , <10" 
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FIG. 9: The angular diameter 9 for the flat model with spin for fi m ,o = 0.3 and f2 Si o = 0.1,0, —0.005 (top, middle, bottom, 
respectively). The minima for these cases are 1.04, 1.605, no minimum, respectively. The spin causes the minimum to move 
right (toward higher z) and the minimum value of decreases. 

Therefore, the term f2 s> o(l + z) 6 , describing the spin effects, is constrained by the BBN because it requires the change 
of expansion rate due to this term to be sufficiently small, so that an acceptable helium-4 abundance is produced. 

VIII. CONCLUSIONS 

We discussed dynamics of the FRW model with spin fluid called the Weyssenhoff fluid. We showed that the effect 
of spin fluid is equivalent to fictitious fluid with equation like the Zeldovich stiff matter and negative energy density. 

The dynamics of these models was analyzed on the two-dimensional phase plane and the Hamiltonian formalism 
was adopted to analyze all evolutional paths in the configuration space. 

The dynamics determines the luminosity distance relationship which is used to test the models using the supernovae 
type la observations. It is shown that the presence of spin fluid matter has no influence on the present acceleration 
rate of the universe, and it is not an alternative to the cosmological constant description of dark energy. 

We showed that the difference between the ACDM model and the model with spin fluid would be detectable if 
the content of spin fluid matter was sufficiently large, e.g., Q Sj o = —0.012. This is possible because supernovae with 
redshift z > 1 should be fainter in the model with spin fluid than in the ACDM model. 

The detailed statistical analysis of supernovae data of sample A (analysis of confidence levels, levels of constant 
X 2 and density distribution of probability) was performed. From this analysis we obtained the limit of the density 
parameter for spin fluid £l Sj0 - At the confidence level of la it is equal to —0.012. 

We also applied the test of the minimum of the angular size of galaxies and showed that it is sensitive to the amount 
of spin fluid matter but difficult to detect. We found that the the presence of spin fluid matter lowers the age of the 
Universe. 



35 



30 



25 



20 



^ „ = 0.004 

s,0 

£1 =0 

s,0 

Q. n = -0.013 

s,0 




FIG. 10: The dependence of the age of the universe to (in units of 10 9 yr) on the parameter f2 mj o for the fiat model with spin 
fluid. 



Moreover we demonstrated that the stronger limit can be obtained from the CMB peak locations using the 
Boomerang data. From the uncertainties of the location of first peak we obtained a small interval for the values 
of spin fluid parameter -1.4 • 10~ 10 < rj s , < -10~ 10 . 

We find the formal analogy between the model with spin fluid and the brane model with dust on a brane. In 
both cases the dynamics equations are formally equivalent. In the early Universe these two kinds of matter scaled 
according the same law. For this reason we also considered the positive values of the spin fluid parameter f2 Sj0 . From 
the statistical analysis we observed that a positive value of spin fluid parameter (brane) is most probable, nevertheless 
a negative value of this parameter is statistically allowed. 

In the near future new high-redshift SN la observations will bring on better data. We expect that they will allow 
us to obtain the estimation of f2 Sj o with a significantly lower error and to restrict the limit for spin fluid. In such a 
case the estimation of r2 Sj0 using SN la data will be worth reconsideration because this method does not depend on 
model assumptions, while the BBN limit on O Sj0 is strongly model dependent. 
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